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The collisionless quasi-equilibrium state realized after the cold collapse of self-gravitating 
systems has two remarkable characters. One of them is the linear temperature-mass (TM) 
relation, which yields a characteristic non-Gaussian velocity distribution. Another is the 
local virial (LV) relation, the virial relation which holds even locally in collisionless systems 
through phase mixing such as cold- collapse. A family of polytropes are examined from a view 
point of these two characters. The LV relation imposes a strong constraint on these models: 
only polytropes with index n ~ 5 with a flat boundary condition at the center are compatible 
with the numerical results, except for the outer region. Using the analytic solutions based 
on the static and spherical Jeans equation, we show that this incompatibility in the outer 
region implies the important effect of anisotropy of velocity dispersion. Furthermore, the 
velocity anisotropy is essential in explaining various numerical results under the condition 
of the local virial relation. 

§1. Introduction 

The astronomical objects in our universe such as galaxies and clusters of galaxies 
are believed to be formed through the gravitational interactions among the mass 
elements composing these objects. These objects originated from tiny fluctuations 
around the almost homogeneous background achieve the virialized state through 
gravitational interactions with the time scale as short as the system's free-fall time. 
However, those with the scale larger than the size of a typical galaxy cannot achieve 
the thermal-equilibrium state, since the time-scale of the equilibrium becomes much 
longer than the age of our universe. These objects are called collisionless, since the 
collisional effect of the gravitational two-body encounters are much less effective than 
that of the gravitational mixing through the potential oscillation. 

Collisionless self-gravitating systems (SGS) eventually settle down to a quasi- 
equilibrium state through the phase mixing and the violent relaxation processes 
under the potential oscillation. 1 ) This quasi-equilibrium state is a prototype of the 
astronomical objects such as galaxies and clusters of galaxies. These objects often 
show universal profiles in various aspects. 

The quasi-equilibrium states for SGSs have been well examined with N-body 
simulations both for the cold collapse simulations starting from the small virial ratio 2 ) 
and for the cosmological simulations starting from the tiny fluctuations around the 
homogeneous expanding background. 3 ) Here we pay our attention to the case with 
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the cold collapse and examine the common characters of the bound state through 
N-body simulations. Since the cold collapse is considered to be a typical case which 
causes the violent relaxation, it is meaningful to examine this typical case in order 
to grasp the universal character of violent relaxation as a first step. 

When the system experiences violent gravitational processes such as cold col- 
lapse and cluster-pair collision, we obtained a universal velocity distribution profile 
expressed as the democratic (=equally weighted) superposition of Gaussian distribu- 
tions of various temperatures(DT distribution hereafter) in which the local temper- 
ature T(r) is defined using the local velocity variance (v 2 ) as T(r) := m(v 2 )(r)/3kB, 
where m is the particle mass. 4 ) Moreover, we have found that the locally defined 
temperature linearly falls down in the intermediate cluster region outside the central 
part, provided it is described against the cumulative mass M r , i.e., dT/dM r = const. 
This fact is consistent with the appearance of DT distributions. 

In addition to the linear TM relation, we have also obtained another peculiar 
fact that the LV relation between the locally defined potential energy and kinetic 
energy holds except for the weak fluctuations, that is, the local temperature T(r) is 
proportional to the local potential $>(r), with constant proportionality 6ksT(r) = 
—m<P(r) for a wide class of cold collapse simulations. 4 )' 6 ) 

In secHJwe will show that the LV relation for SGS is supported by the results 
of a variety of cold collapse simulations and that among the spherical and isotropic 
models, polytropes with index n ~ 5 with a flat boundary condition at the center 
are compatible with the numerical results except for the outer region of the bound 
state. In secEl we examine the anisotropic models under the condition of the LV 
relation. We will show that the analytical solutions exist under the LV condition 
with constant (3. These analytical solutions can be utilized to form the simple model 
with the anisotropic velocity dispersion. Sec|I] is devoted to the summaries and 
conclusions of this paper. 

§2. LV relation and analysis with isotropic models 

In the previous paper, we showed that linear TM relation is one of the remarkable 
characters of SGS bound states after a cold collapse. 4 ) Here we pay our attention to 
another characteristic for SGS bound states, that is, the LV relation. 

It is well known that the gravitationally bound system approaches a virialized 
state satisfying the condition 



where W and K are, respectively, the averaged potential energy and kinetic energy 
of the whole bound system. This is a global relation that holds for the entire system, 
after the initial coherent motion fades out. 

Here we define the locally averaged potential energy and kinetic energy inside 
the radius r, respectively as 



W + IK = 






Local virial relation and velocity anisotropy for collisionless SGS 



3 



K r 



(v 2 )(r>) 



p(r')47rr' dr', 



(2-3) 



where *(r') means the local object * evaluated at r' . 
Then we extend the virial relation (|2-1() locally as 

W r + 2K r = 0, (2-4) 

and examine how precisely this relation is locally attained inside a bound state. 



1.5 



-a 
0.5 



I H 

SC(b,=0,a=0) 

SC ( b t =0.1, a=0 ) hBh 

SC( b t =0.3,a=0) i — Q — i 

SC(b t =0.5,a=0) mAh 



SC( b t =0, a=l) 
SC( b t =0, a=2) 
CC 



0.2 0.4 0.6 0.8 1 

M r /M tof 

Fig. 1. The LV relation for some numerical simulations obtained by a typical cold collapse simu- 
lation; spherical collapse (SC) and cluster-pair collisions (CC). In the case of SC, 5000 particles 
are distributed with a power law density profile (p oc r~ a ) within a sphere of radius R and the 
initial virial ratio (bt) is set to be small. In the case of CC, each cluster has the equal number 
of particles (2500) and all particles are homogeneously distributed within a sphere of radius 7? 
and is set to be virialized initially. The initial separation of the pair is 67? along the x axis. In 
all of the simulations, softening length e = 2~ 8 R is introduced to reduce the numerical error 
caused by close encounters. The local virial ratio b is plotted as a function of M r /Mtot, where 
Mtot is the total mass of the system.. The virial ratios at M r are time averaged from t = 5tff 
and t = WOtff. 

More precisely, the above relation Eg. 1)2-4(1 is equivalent to the purely local re- 
lation 

2(v 2 )(r) = -#(r), (2-5) 
at each position r. Hence, we define the LV ratio 

b(r) = -2(v 2 )(r)/$(r), (2-6) 

and examine the value of b(r) for each shell. For a wide class of collapses including 
cluster-pair collision, 4 ) the value b(r) takes almost unity: it deviates from unity less 
than 10 percent upward in the intermediate region and a little more downward in 
the inner and outer region for all of the simulations (Fig^). Hence a wide class of 
cold collapse simulations yield the relation (J2-5() quite well. This LV relation should 
be another characteristic of SGS, in addition to the linear TM relation. 
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From a viewpoint of these two characteristics, we compared our results with 
polytropes. The polytrope with n = 5 is special among those models, since it 
exactly satisfies the LV relation (Fig|2J). Moreover, it admits the analytical solution 
called Plummer's model under the condition that the potential is flat at the center 
of the system. 5 ^ We got the result that Plummer's model satisfies the linear TM 
relation quite well in the intermediate region 0.2M tot < M r < 0.8M to t- 6) On the 
other hand, both its TM relation and density profile in the outer region differ from 
the ones derived in cold collapse simulations: the temperature falls off more steeply 
and the density behaves as r -5 in the outer region. 
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Fig. 2. LV relation for polytropes. The solid curve refers to the LV ratio b(M r / M to t) as a function 
of cumulative mass M r /Mtot for Plummer's model. The other curves refer to b(M r /M to t) for 
polytropes with n = 3.0 (dashed), n = 1.5 (dot-dashed), and n = 0.51 (dotted). 

§3. Models with anisotropic velocity dispersion 

In the previous section, we commented that Plummer's model expresses the 
bound state quite well after a cold collapse among the spherical systems with isotropic 
velocity dispersion. However, it fails in explaining the outer region of the bound state, 
where the velocity dispersion is anisotropic (FigEJ). Here we pay our attention to this 
anisotropy and examine the general static Jeans equation admitting this anisotropy. 

In the collisionless system, the particle distribution can be described by the 
Vlasov equation. Integrating the spherical and static Vlasov equation in velocity 
space after multiplying v T results in the spherical and static Jeans equation, 

d{p(v 2 r )) # 2(v 2 r )PP 

ar ar r 

where (3 = 1— ((v 2 ) + (t^)) / {^\ v r )) ^ s the anisotropy parameter and (p is the relative 
potential defined as <p = ~ ^> where ^* is the maximum energy of the particles 
in the system. 5 ) 

Assuming that the density is non-zero everywhere but the total mass is finite, 
we can get the condition that <P* = 0. In this case, from the condition that (v 2 ^ = 
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Fig. 3. The distribution of the anisotropy parameter [3(r). A plot with error-bar is the result of 
the numerical simulation SC with (bt,a) — (0,0) but with TV = 10000 in FigQ A solid thick 
line represents the combined solution (IA- It with c — 0.2 and a dot line represents the analytical 
solution with (3 — 0.5. 



(3 - 2(3) (v 2 ), the LV relation 1(73)) takes the form of 



~<t>, 



where r = 1/(3 — 2(3) . From Eqs. (|3-1|) and ()3-2|) . we derive the relation 

d (rpcj)) _ _ (3r - 1) p4> 

dr dr r 

The relative potential also follows the Poisson's equation 

]_d_ f 2 df \ = _ 47T 2 Gp 
r 2 dr \ dr 1 



(3-2) 



(3-3) 



(3-4) 



Eas. (|3-3j) and ()3-4(l are basic equations for cj) and p when the anisotropy parameter 
j3 is given. 

The equation (|3-3[) is integrable under the condition that [3 is constant, which 
leads to the relation 

p = A^-^r- 2 ! 3 , (3-5) 

where A is an integral constant. Substituting Ea. H3-5|) into Ea. (|3-4j) . we can get the 
nonlinear equation for (f>, which admits the analytical solution 
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(3-6) 



where s = 2(1 — 0) and ro = (^f^p) M 2 ot under the boundary condition at infin- 
ity; 



<p (oo) = p (oo) = 

_ _GM to t 

dr 



(3-7) 
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. 7 ) We will analyze more general solutions for eq. (|3-1[) and 1)3-4(1 under the constant 
LV ratio b in another paper. 8 ) 

Especially when (5 = 0, the solution ()3-6() is nothing but Plummer's model. The 
Plummer's model explains several characters of the bound state after a cold collapse 
or cluster-pair collisions quite well except for the outer region. Here we compare 
this solution with the solution connecting inner Plummer's solution with the outer 
analytical solution with constant (3 (See AppendixEJ. Of course, the actual function 
form of (3(r) is too complicated to be expressed by the connection of several constant 
(3 analytical solutions. However, we can roughly see that considering the anisotropy 
reduce the incompatibility of the solution to the bound state after a cold collapse 
(FigHand Fig©. 
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Fig. 4. A log- log plot for a density profile with the unit of Vh = Mtot = G = 1, where th is the half- 
mass radius of the bound system. A plot with error-bar is the result of the numerical simulation 
SC with (bt,a) = (0,0) but with N = 10000 in Fig0 A solid line represents the combined 
solution lA-lll with c = 0.2. A dot and dot-dashed line represent the analytical solution with 
/3 = 0.5 and the Plummer's model, respectively. The coincidence between numerical data and 
the combined solution is remarkable. 



§4. Conclusions 

The bound states after a cold collapse or several other initial conditions with 
small initial virial ratio are characterized by the linear temperature-mass (TM) rela- 
tion and the local virial (LV) conditions. Among the spherical and static models as 
the solution of Vlasov equation with isotropic velocity dispersion, the polytrope with 
n = 5 is very special, since it satisfies the LV relation exactly. Plummer's solution, 
which is the analytical solution of polytrope with n = 5 satisfying the flat boundary 
condition at the center explains several characters of the bound state quite well ex- 
cept for the outer region. Here we showed that this incompatibility can be reduced 
by adding the anisotropic effect under the condition of the LV relation. This leads 
to the result that the bound state after a cold collapse characterized by the density 
profile with the asymptotic form oc r -4 and linear TM relation can be explained 
quite well by the anisotropic model satisfying the LV relations. 

Here we have restricted our analysis to cold collapse models with small initial 
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Fig. 5. A Temperature-Mass relation with the unit of rj» = Mtot = G = 1, where ru is the half- 
mass radius of the bound system. A plot with error-bar is the result of the numerical simulation 
SC with (bt,a) = (0,0) but with N — 10000 in Fig^ A solid line represents the combined 
solution lA-lll with c = 0.2. A dot and dot-dashed line represent the analytical solution with 
P = 0.5 and the Plummer's model, respectively. The coincidence between numerical data and 
the combined solution is remarkable. 

virial conditions. However, it seems important how widely our results are applicable 
in more general mixing processes including merging process in cosmological simula- 
tions. It is well known that the bound states observed in cosmological simulations are 
quite different from those in cold collapse. For example, the density has the central 
cusp oc r~ x and takes the outer asymptotic form oc r~ 3 in cosmological simulations. 3 ^ 
Moreover, the phase space density defined as p/a 3 has a single power-law property 
in radius. 9 ) Hence it seems important to see if the LV relation is admitted even in 
the case with such different characters. These points are now under investigation. 8 ) 



Here we connect the two analytical solution at r = r c ., i.e., Plummer's model 
((3 = 0) for r < r c and (3 = 0.5 analytical solution for r > r c . Defining <p c = <f> (r c ) and 
describing the physical quantities with the unit of r c = (j) c = G = 1, the connected 
solution is described with one-parameter c as 



Appendix 





from the continuity condition of (p(r) and d<j)(r)/dr at r=l 





(A-2) 



Hence the total mass M to t is described as 



M tot 



lim M (r) 



1 +3c 2 . 
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The mass ratio between the total mass of inner Plummer region, M p i um and Mtot 
becomes 

Mpi um ^ 1 
Mtot ~ (l + 3c 2 ) 2 ' 

Hence the parameter c determines the mass ratio of the connected two regions. 
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